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We consider magnetohydrodynamical interaction between relativistic pulsar wind 
^ | and static magnetosphere in binary pulsar system PSR J0737-3039. We construct semi- 

analytical model describing the form of the interface separating the two pulsars. An 
' assumption of vacuum dipole spin down for Pulsar B leads to eclipse duration ten times 

■ longer than observed. We discuss possible Pulsar B torque modification and magnetic 

field estimates due to the interaction with Pulsar A wind. Unless the orbital inclination 
. is < 86° , the duration of eclipses is typically shorter than the one implied by the size of 

the eclipsing region. We propose that eclipses occur due to synchrotron absorption by 
t^J- | mildly relativistic particles in the shocked Pulsar A wind. The corresponding optical 

depth may be high enough if Pulsar A wind density is at the upper allowed limits. 
We derive jump conditions at oblique, relativistic, magnetohydrodynamical shocks and 
discuss the structure of the shocked Pulsar A wind. Finally, we speculate on a possible 

8 ■ mechanism of orbital modulation of Pulsar B radio emission. 
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1. Introduction 

A recent discovery of eclipsing binary pulsar system PSR J0737-3039 (Burgay et al. 2003; 
Lyne et al. 2004) may serve as an important tool in studying close environment of the neutron 
stars. In this system a fast recycled Pa = 22.7 msec Pulsar A orbits a Pb = 2.77 sec Pulsar B on a 
2.4 hour orbit inclined at ~ 87.7° (with uncertanity of the order of unity) to the line of sight. Pulsar 
A shows 27 sec frequency independent eclipses at the moments of superior conjunction (Kaspi et 
al. 2004). In addition, Pulsar B shows orbital-phase dependent variations of intensity. It is also 
very weak at the inferior conjunction, consistent with being eclipsed as well (S. Ransom, private 
communication) . 

In this paper we first construct a geometrical model of Pulsar A eclipses. Similarly to Arons 
et al. (private communication) we propose that the interaction of Pulsar A wind with Pulsar B 
magnetosphere resembles Solar wind - Earth magnetosphere interaction (see Fig 1). As a result, a 



o 
S3 

Oh 



- 2 - 



bow shock forms in Pulsar A wind, separated from Pulsar B magnetosphere by a sheath of shocked 
material. 

Absorption of Pulsar A radio emission may occur either in the shocked Pulsar A wind or in 
the magnetosphere of Pulsar B. If absorption occurs at cyclotron resonance with non-relativistic 
particles present on the closed field lines, then it is expected to be strongly frequency dependent, 
which is not observed. Synchrotron absorption by relativistic particles on the open field lines can 
give a frequency independent eclipse. Still, we disfavor this possibility, First, there are no indication 
of modulation of the eclipse by the rotation of Pulsar B (though the present temporal resolution 
may not be sufficient). Secondly, since the line of sight passes at the edge of magnetosphere (see 
Fig. 3), and since we see pulses from Pulsar B, it is hard to imagine geometry in which near the 
superior conjunction the line of light constantly passes through open field lines, giving full eclipses. 
Thirdly, the dimming of Pulsar B near the inferior conjunction is consistent with absorption being 
due to the shocked Pulsar A wind. 

An alternative possibility for the location of the absorbing material, which we favor, is that 
the absorption occurs in the shocked Pulsar A wind, enveloping Pulsar B. Though the size of the 
shocked sheath is expected to vary by about 25% due to Pulsar B rotation, it may still have an 
approximately constant column density and optical depth regardless of its size. 

What determines the size of the eclipsing region? It may be determined either by geometric 
factors, so that the beginning and the end of the eclipse occur at the moment when the line of sight 
enters the absorbing region, or by microphysical processes in the plasma so that the eclipse occurs 
at the moment when optical depth r to scattering/absorption becomes of the order of unity. We 
strongly prefer the geometrical interpretation of eclipses since any relevant microphysical process 
gives a frequency dependent absorption coefficient so that the condition r ~ 1 would occur at 
different points for different frequencies. 

The geometrical (as opposed to microphysical) interpretation of the eclipse duration seems, 
at first sight, inconsistent with the data: a simple scaling of the size of Pulsar B magnetosphere 
would predict an eclipse which is much longer than observed. The light cylinder radius of Pulsar 
B is 1.29 x 10 10 cm, so that its angular size seen from Pulsar A is 0.14 rad. The line of sight passes 
2.3° = 0.04 rad = 3.6 x 10 9 cm away from Pulsar B. Eclipsing part, corresponding to a 27 sec 
eclipse duration, is 1.8 x 10 9 cm Kaspi et al. (2004), 1 which subtends a total angle of 0.02 radians 
(see Fig. 2). Thus, the size of the eclipsing region is ~ 0.04 rad = 4 x 10 9 cm, which is more than 
3 times smaller than the light cylinder radius of Pulsar B. 

A possible resolution of this disagreement is that the pressure of Pulsar A wind compresses 
Pulsar B magnetosphere to lateral sizes much smaller than the light cylinder. In this paper we 
first calculate semi-analytically the form of the interface between Pulsar A wind and Pulsar B 



x At this point magnetic field of the B pulsar Bb = 13 G, so that the cyclotron frequency vb = 4.4 x 10 7 Hz, which 
is too low to produce cyclotron absorption if particles are non-relativistic. 
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magnetosphere. We find that if the magnetic field of Pulsar B is calculated using the vacuum dipole 
formula, the lateral size of the sheath is still much larger than the one inferred from eclipse duration. 
On the other hand, electromagnetic interaction of Pulsar B with the sheath strongly modifies the 
structure of the magnetosphere, so that the vacuum dipole formula becomes inapplicable. The 
modification of the spin down torque leads to an estimate of Pulsar B magnetic field which is 
typically 3-5 times smaller. Still, the size of the sheath is much larger than inferred from eclipse 
duration. A possible resolution is that inclination angle is < 86°, so that the line of sight cuts at 
the edge of the sheath. 

As independent parameters of the system we chose the spin down luminosities La = 5.8 x 10 33 
erg/sec and Lb = 1.6 x 10 30 erg/sec and the separation of D = 9 x 10 10 cm. These are the most 
reliable parameters, while the quantities like the inferred surface magnetic fields depend on the 
details of the magnetospheric structure. 



2. Model of Pulsar A wind — Pulsar B magnetosphere interaction 

Consider a point C on the interface between Pulsar A wind and Pulsar B magnetosphere 
defined by the radius- vector tb(0b) (see Fig. 2). A normal to the interface makes an angle 

tana = <% s mr^^e) (1) 

with the radius- vector tb{0b)- Here 9b is the polar angle between r# and the normal to the plane 
of the orbit. At the same point, the radial flow from Pulsar A, 2 leaving at the polar angle 6a, is 
inclined at angle 4> to the interface. Angles a, 4>, 9a and 9b are related by 4> = 7r /2 + a — (9a + Ob)- 
From the triangle ABC we also find 

tan6U + tan6> B = (2) 

tb cos Ob 

and 

r-A = cos Ob 
vb cosOa 

We assume that the position of the interface is determined by the pressure balance between 
the static pressure of Pulsar B magnetic field and the dynamic pressure of Pulsar A wind (this is 
often called Newtonian approximation) : 

L A sin 2 <p _ ji 2 



2 Pulsar A wind may be assumed radial since the wind velocity, which is of the order of the speed of light, is 
much larger than the orbital velocity, ~ 300kms _1 and much larger than the angular rotational velocity of the wind, 
v ~ c 2 /(DQa) ~ 0.001c, where Qa = 2-k/Pa is the angular rotation frequency of Pulsar A. 
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where \x is the magnetic moment of the B pulsar and /(Ob) depends on the inclination of the 
magnetic moment. /(Ob) = 1, 1 + Ss'm 2 Ob, 1 + 3cos 2 #b for magnetic moments oriented along 
x, y and z correspondingly. 

Next, we dimentionalize the problem by introducing r# — > rD, £ = 2LaD 4 /^ 2 , 

/(Ob) = £ r r + c0Bd B do B lnr B (0B) -smO B 

r 6 l + (de B \nr B (0 B )) 2 \ l-2rsinfl B + r 2 

which for a given /(Ob) is an ODE for t(0b) depending on £ as a parameter. Next, instead of £ we 
introduce a dimensionless (measured in terms of D) stand-off distance tq 

f = ^-^/(*/2) (6) 

and parametrize r — > rog(0). This gives 

/(gg) _ ff 6 / (l-ro)(r g 2 -gsin^ + ff / cosg B ) 
/(tt/2) 5 2 + 5 ' 2 V 1 + r 2 <7 2 - 2r o5 sin B 

For a given stand-off distance ro Eq. (7) determines the form of the interface. We integrate this 
equation numerically, limiting integration to the regions close to Pulsar B light cylinder (near the 
apex point g ~ 1 + 2(7r/2 — 0b) 2 /3). Further down the stream the approximation of a static 
magnetosphere of Pulsar B breaks down and the balance at the contact discontinuity is determined 
by the pressure balance of two winds. 





3. Magnetic field of Pulsar B 

In order to find the stand-off distance we need to estimate a magnetic field strength of Pulsar 
B. This is not straightforward since the conventionally used vacuum dipole formula is not applicable 
in this case, as we argue below. 



3.1. Vacuum dipole spindown 

Conventionally, magnetic field of pulsars is determined by the vacuum dipole spin down formula 



L = |< s < sC (^) 4 (8) 

Even in case of isolated pulsars this is only an approximation since a large contribution to the 
torque comes from the currents flowing on the open field lines of magnetosphere. ( If the typical 
current density is the Goldreich- Julian density times velocity of light, then the torque from the j x B 
force integrated over the open field lines gives the same estimate as the vacuum dipole luminosity 

(8)0 
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Equating (8) to the spin down luminosity Lb = fltll, where / ~ 10 45 gcm 2 is the moment of 
inertial of neutron star, gives Bb = 1.5 x 10 12 G (Lyne et al. 2004). In this case £ = 4 x 10 7 , 
vq = 8.6 x 10 -3 for /(tt/2) = 1 (when the B pulsar magnetic moment is oriented along x or z axis), 
and ro = 1.07 x 10~ 2 for f(ir/2) = 4 (when the B pulsar magnetic moment is oriented along y axis). 
The corresponding stand-off distances are r m = 4.8 x 10 9 cm and r m = 6 x 10 9 cm. 

For these parameters ro we integrate Eq. (7). The results are given by upper set of curves in 
Fig 3. It is evident that the resulting size of the interaction sheath is much larger than the observed 
duration of the eclipse (the expected duration of the eclipse is ~ 10 times longer than observed). 



3.2. Modification of the spindown of Pulsar B due to interaction with Pulsar A 

wind 

Applicability of the vacuum spin down formula (8) in the case of interacting system PSR J0737- 
3039B is even more doubtful than in the case of isolated pulsars. Since the magnetospheric radius is 
located deep inside the light cylinder, the structure of Pulsar B magnetosphere is strongly distorted 
if compared with the vacuum case. The details of the structure are bound to be complicated. 

Qualitatively, there are three possibilities for the structure of Pulsar B magnetosphere and 
corresponding spindown torque, (i) the interface may be a perfect conductor so that no Pulsar B 
field lines penetrate it. (ii) the interface may be a perfect resistor, so that all the Poynting flux 
from Pulsar B reaching it gets dissipated, (iii) the interface may be partially resistive so that large 
surface currents are generated, which together with the poloidal field of Pulsar B produce a torque 
on the star. 

Modification of the spindown torque in all these cases are very different. We are mostly 
interested in cases (ii) and (iii). For case (ii), assuming that the typical current on the open field 
lines is of the order of the Goldreich-Julian current and that the size of the open field lines is 
determined by the magnetospheric radius r m , the spin down becomes 

2 / \ 2 

(9) 



t = (1)^(0^) (^ 



Equating (11) to spin down luminosity Lb and using the force balance (4), we find 

3 3 / 4 VOLT „ 
Bb = -= tja ^— = 4.7 x 10 11 G 




V2\l aJ Vn < 10 > 

Using parameters (10) we integrate Eq. (7) (the middle solid curve on Fig. 3; in this case ro = 0.045 
for /j, = fi x ). The size of the sheath is only slightly modified if compared with the vacuum dipolc 
and is much larger than the one implied by eclipse duration. 
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Finally, the interface may be partially resistive. Relativistic boundary layers are strongly 
unstable on scales of tens to hundreds of gyroradii (Smolsky & Usov 1996; Liang et al. 2003), so 
that for ultra-relativistic electrons with 7 ~ 10 6 the kinetic thickness of the interface may become 
macroscopically large. As a result, Pulsar B field may penetrate the interface, similar to what may 
happen in neutron star - disk interaction (e.g. Ghosh &Lamb 1979a,b; Wang 1995). In this case 
the poloidal magnetic field of Pulsar B will be twisted by the pulsar rotation to produce toroidal 
magnetic field which may be as high as poloidal magnetic field at the interface, B^(r m ) < BB(r m ). 
(Since the magnetospheric radius is much smaller than the light cylinder radius, light travel effects 
may be neglected.) The resulting surface current on the interface together with the poloidal field 
will produce a j x B force and a torque on the pulsar. The spin down then becomes 

A factor 1/2 has been introduced in front to account for the fact that only the part of pulsar B 
magnetosphere facing Pulsar A experiences the torque. 

Equating (11) to spin down luminosity Lb and using the force balance (4), we find 
B B = A = 3 x 10 11 G 

3\ 1/3 (L B V' Z fcD^ 1/3 



L=[^)[- )B^ NS c[^\ l^f.) (11) 



2j \L A J V ^ 
Magnetic field in this case is ~ 5 times smaller. 



3.5 x 10 9 cm (12) 



Using parameters (12) we integrate Eq. (7) (the lower solid curve on Fig. 3; in this case 
ro = 0.038 for fx = fx x ). 

Both modified magnetic field estimates (10-12) predict the size of pulsar B magnetosphere 
which is 2-3 times larger than inferred by eclipse duration. A possible resultion of the disagreement 
is that the orbital inclination angle is < 86°. 



3.3. Spin down due to propeller effect and Magnus force on Pulsar B 

Interaction of Pulsar B with the wind of Pulsar A may also produce propeller effect, whereby 
the material of Pulsar A wind is flung off by the rotation of Pulsar B. As a result, there will be 
extra torque on Pulsar B produced by the Magnus force (a force due to a difference in pressures 
at the two sides of the Pulsar B magnetosphere with rotation velocity aligned and counter-aligned 
with the radial wind of Pulsar A, c.f. a "dry leaf kick in soccer). Next we estimate the importance 
of the propeller effect on the spin down of Pulsar B. 

The torque of the Magnus force acting on Pulsar B can be estimated as 

1 ~ 4W [L6) 
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Equating this to the change of angular momentum we find the corresponding magnetospheric radius 

r ™ - 22/5 (£) 1/5 " 5 = 7 ^ 109 - (") 

This radius is larger than the one given by the dipole formula. Thus, we conclude that the propeller 
effect is not important for the spin evolution of Pulsar B. 3 



4. Eclipsing mechanism: synchrotron absorption 

We propose that the absorption of Pulsar A radio beam occurs due to the synchrotron absorp- 
tion in the shocked plasma. We parametrize Pulsar A wind by a magnetization parameter a, the 
ratio of Poynting to particle fluxes (Kennel & Coroniti 1984a). The magnetic field in the sheath 
(in the laboratory frame) is 

b =VtI^I =2 Vt^ g=1 - 2G <15> 

where the last equality assumes a = 3 x 1CT 3 , a value inferred for Crab pulsar (Kennel & Coroniti 
1984a). The factor of 3 in front assumes a compression in strong relativistic fluid shock (see 
also Section 5 for more details). The non-relativistic cyclotron frequency inside the sheath is 
vb = 2.4 MHz. In order to absorb at the observed frequencies u f) S ~ 1 GHz, the particles should 
be relativistic with a Lorentz factor 7 ~ \J u b s /uB = 30. It is not obvious at all that such particles 
are present in the shocked flow. Estimates of the bulk Lorentz factor in case of the Crab pulsar 
give T ~ 10 6 (Kennel & Coroniti 1984b). If a similarly high Lorentz factor is assumed for PSR 
J0737-3039A, then one expects that the lowest energy cut-off in the shocked flow is similarly high, 

7 ~ r. 

On the other hand, radio emission of the Crab nebula is attributed to electrons with much lower 
energies, 7 <C T. These electrons may be either cooled remnants of the intense early injection Atoyan 
(1999), or, more likely, may represent a completely different population of accelerated particles. 
Thus, invoking the Crab pulsar as an example, it is feasible that relativistic magnetized shocks do 
produce a low energy population of electrons. Note, though, that since the magnetic fields in the 
bulk of the Crab nebular and in the sheath of the interacting winds of PSR J0737-3039 differ by 
some 4 orders of magnitude, the corresponding Lorentz factors of synchrotron emitting/absorbing 
particles differ by two orders of magnitude. Unfortunately, in the absence of understanding of the 
acceleration mechanism of radio electrons we cannot judge whether this is possible. (Alternatively, 
it is feasible that the low energy population in the sheath comes from Pulsar B due to a "leakage" 
through the interface.) 



3 The same Magnus force also produces a torque on the orbital motion of Pulsar B. The corresponding effect on 
the orbital evolution is too small to be of any importance. 
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Below we assume that a population of low energy electrons with a power law distribution n(e) = 
K€~ p is indeed present in the sheath (e is the energy of particles). The parameter k is related to the 
density of pairs and the low and high energy cut-offs in the spectrum, k ~ nmc 2 / (p— l)7^~n ~ nmc2 
where we assumed that ^min ~ 1 an d p close to 2. 

Cyclotron absorption coefficient is then (Lang 1974) 

q° = 2 x 1(T 2 g{p) (3.5 x 10 9 ) p k B ip+2 ^ 2 v'^)/ 2 cm" 1 (16) 

where g{p) is a coefficient of the order of unity, v is observing frequency. We also assumed that 
radiation propagates orthogonally to the field lines. 

We normalize the pair density in the wind to the Goldreich-Julian density at Pulsar A light 
cylinder. The particle flux for the Pulsar A is then Na = X^AttLac/c 2 ~ A x 10 32 s _1 and the 
density at distance of Pulsar B is ua = NaU^cD 2 ) = 3 x 10~ 2 A cc. In order to estimate the 
thickness and the particle density in the sheath it is necessary to know the details of the flow 
structure. As a simple estimate we assume that the column density through the sheath is of the 
order of the column density of Pulsar A at the distance of Pulsar B, ~ uaD. For particle index 
p = 2 we find 

Thus, in order to produce optical depth > 1 the multiplicity factor should not be quite large 
A > 10 6 (flatter spectra p < 2 give higher optical depth also, for p = 1.5 it increases by an order 
of magnitude). The required multiplicity is large, but not impossible (a value of multiplicity factor 
invoked for the Crab is A = 10 6 , (Kennel & Coroniti 1984a), see also Arons & Scharlemann (1979); 
Muslimov k Harding (2003)). 

We conclude that synchrotron absorption by weakly relativistic particles in the shocked Pulsar 
A wind is a possible eclipsing mechanism, but acknowledge that the required wind multiplicity is 
at the upper possible end. 

5. Dynamics of the shocked Pulsar A wind 

In the previous section we calculated semi-analytically the form of the interface between Pulsar 
A wind and Pulsar B magnetosphere. In reality, the interface will consist of the bow shock and a 
contact discontinuity (see Fig 1). 4 In this section we study the dynamics of the shocked pulsar wind 
in the sheath. We first derive the oblique jump conditions for relativistic MHD shocks. Relativistic 
oblique fluid shocks have been considered by Konigl (1980); we are not aware of a work which 
considered relativistic oblique MHD shocks. We apply the results to PSR J0737-3039, assuming 
that the form of the shock follows the form of the interface calculated in Sections 2 and 3. 



4 In fact, since magnetic field of the Pulsar A wind piles up on the contact discontinuity, the two media are separated 
by a rotation ( Alfven ) discontinuity. 
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5.1. Oblique fluid shocks 

Let the stream lines make an angle <f> with the shock and let the post-shock flow make an angle 
X with the initial velocity (Fig. 5). Oblique shock conditions can be obtained from normal shocks 
and an additional condition that the component along the shock velocity remains constant. 

n\U\ sin <f> = ri2U2 sin(x — (f>) 

wiu\ sin 2 (f) + pi = w 2 ul sin 2 (x - <j>) + P2 

wijiui sm<f> = W2J2U2 sin(x - 4>) 

uj cos 2 4> = u 2 2 cos 2 (0- x) , lg x 
1 + uf sin 2 1 + u 2 sin 2 (<f> - x) 

where w is enthalpy, u is four-velocity, n is density, p- pressure and 7 is Lorentz factor; velocity of 
light has been set to unity. Subscripts denote unshocked (1) and shocked (2) fluids. Relations (18) 
can be resolved (Landau & Lifshitz 1975) 



u\ sin i 



(e 2 +Pi)(V2 ~Pi 



(ei +Pi)(e 2 - ei - (p 2 - Pi) 



u 2 sin(x - 4>) 



(ei +P2){p\ -P2) 



(e 2 + p 2 )(ei - e 2 - (pi - P2) 



ni _ u 2 sin(x - 4>) _ (ei +j? 2 )(ei + Pi) _ ^ ^ 
n 2 uisincp y (e 2 + Pi)(e 2 + P2) 

where w = e + p and 77 is the compression ratio, defined here as a ratio of rest frame densities. 

For adiabatic flow, using relations 

^ = n+-^-p=^-pr a (20) 

where u s is a sound four-velocity, and T a is adiabatic index. The compression ratio 77 can be 
expressed as a function of u S: \ and u S)2 . The corresponding relations are too bulky to be reproduced 
here. In the case of initially cold plasma (pi = 0, w± = e\ = n\, u Sj i = 0) we find 

, (r a -i)(r (i-7/)-(i + 7/)) 



r a (2 - r o (l - V )) V 2 sin 2 
sin 6 



u 2 = TJUl 



sin(x 



„2 (r a -i)(r a (i-r/)-(i + >7)) 
s ' 2 r a (3- 17) -1^(1-1/) -2-1/ 1 j 

In the non-relativistic limit b<1 these relations give 77 = (r a — l)/(r a + 1) + 0(u 2 ), while in the 
strongly relativistic limit u\ S> 1, 77 <C 1, 



77 



r -i 1 



ui sin0y / r a (2 - r a ) 2\/2«isin0 
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r a — 1 1 sin (J) 

y/T a (2 - r a ) ~ 2V2sm( X -4>) 

^-S^.-ji (22> 

where the second equalities assume T a = 4/3 (this is also assumed for all numerical estimates 
below). Thus, for cold flow < r) < (T a - l)/(T a + 1). 

As independent parameters of the problem we chose the initial four-velocity u\ and the com- 
pression ratio rj. For a particular case of cold plasma u\ and rj are related by Eq. (21). 



Eliminating u 2 from Eqns. (18), we can determine x(<fi, u i)'- 

1 + u\ sin 2 (j) 
1 + rp"u\ sin 2 <f> 



9 . , , 9 1 + u 2 sin 2 6 9 , 
tan 2 (0 - x) = V 2 — * . 2 \ tan 2 <j> (23) 



which, using Eq. (21) gives 

T a (l - rj) tan 
l + (l-r a (l-r ? ))tan 

For strong shocks, r/ — > 0, 

T a tan 
1 - (r a - 1) tan 

For T a = 4/3, the maximum deflection angle Xmax = k/6 is reached at 4> ma x = tt/3 



We can also determine the post-shock Mach number: 



> u 2 . u? cos 2 + (1 + nf )r? 2 sin 2 </> 2 cos 2 + 1/4 
M 2 - ~T~ ~ 771 T~, — 2 ■ 2 2 ~ ~ 2~7 V^oJ 



t2 _ 

S u| s 1 + sin^ sin^ 



For initially cold flow, the compression ratio r] can be eliminated using relations (21); then Eq. (26) 
can be solved for the angle <\>m at which the post shock flow is sonic. We find (f)M = <t>max = vr/3. 



5.2. Oblique MHD shocks 

Next we find jump conditions for strong, ultra-relativistic fast magnetosonic shocks, assuming 
that the magnetic field lies in the plane of the shock. (When the field is in the plane of the shock, 
the relevant MHD expressions can be obtainable by generalizing the hydrodynamical relations by 
substituting for the pressure p and internal energy density e: p — > p + b 2 e — > e + 6 2 /87r.) The 
shock jump conditions for relativistic MHD shocks are Landau & Lifshitz (1975) 

biUi sin<j) = b 2 u 2 sin(x — (ft) 

b 2 

(wi + bj) u\ sin 2 d> + Pl + -± = (w 2 + bf) u\ sin 2 (x 
{w\ + bf) u\ sin 07i = (w 2 + bf) u 2 sin(x - <f>)^ 2 



+P2 + i 



(27) 
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(continuity equation and equation for velocity along the shock remain the same). In Eq. (27) b is 
a rest frame magnetic field times y/47r. Equations (27) can be resolved 



Ui 



((&? - bf)/2 + pi - p 2 ) {{b\ + bl)/2 + pl-p 2 + w 2 ) 1 
2{b\ +u/l)(2(p 2 -Pi) + w2 - wl) sin0 



sm <p 
sin(0 - x) 



= ni = h = n 2 sin((/»- X ) = / (^i + frf) (wi - pi + p 2 + (&j + gVg) 
77 ra 2 & 2 sin V (u>2 + & 2 ) (™2 - P2 + Pi + (&? + &!)/ 2 ) 

To characterize the magnetization of the flow we introduce parameter (see Kennel Sz Coroniti 
(1984a)) 

* = (29) 

Wl 

which is the ratio of the rest frame magnetic and plasma energy densities (it is also equal to the 
ratio of Poynting to particle fluxes). 

Similarly to the fluid case, the compression ratio (28) can be written in terms of preshock and 
post-shock fast magnetosonic four-velocities 

2 pT a + wa 

u„ = — (30) 

If the preshock plasma is cold and in the limit of strong shocks rj — > (it is also necessary that 
u\ ^> \f&)-, and assuming T a = 4/3, we find 

= 1 2 + T a (a -2) -4a- y/Tl(2 + a) 2 + 4(1 + 2a) 2 - 4r a (2 + 3a + 2a 2 ) 

Ul ~ * \| 8(2 - r a ) ( y/Tl(2 + a) 2 + 4(1 + 2a) 2 - 4r a (2 + 3a + 2a 2 ) + 2 + (4 - T a )a) " 

1 + 4a + VI + 16a(a + 1) _ f(a) 



2V2 V ^3 + 4a + y/1 + 16a(a + = fj sin< ^ 



2 - r a v 2 



^-JTE^E-Jl^. (31) 



(cf. Kennel & Coroniti (1984a)). Equation (31) defines the function /(a). 

Equations describing oblique relativistic MHD shocks remain the same as in the fluid case, 
with dependence u\(rj) given by (31). The deflection angle becomes 



(/M + ^/TTTEflw (32) 

v / TT7W I -/Mtan 2 * 
The maximum deflection angle is reached 

cos 4> max = ^ f(a)(^/TTJW - f(^)) (33) 
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As a -> oo, </> ma:E -> vr/4, Xm«x -> V(4cr). 

Similarly to the hydrodynamic case, using Eq. (26), we can find the angle </>i at which the 
post shocked flow remains sonic. In the ultra-relativistic limit we find 

. / (2-r Q )(l + /(q)^ 

sin <t) M = y (34) 

(see Fig. 7). As cr — > oo, 0m — ► arcsin = 54.73°, xm 1/(3\/2<t). 

Using the relations derived in the previous section we can calculate the post shock velocities 
(Fig. 8). At larger angles the post-shock flow becomes relativistic, T s > 1 (T s is the post-shock 
Lorentz factor). 



6. Pulsar B emission 

Another puzzling property of PSR J0737-3039 is the variations of Pulsar B flux depending 
on the orbital phase (Lyne et al. 2004). In addition to two sections of the orbit where it is very 
bright, so that single pulses can be seen, the flux density of the pulsar B emission is at a minimum 
near, but centered slightly before, the inferior conjunction. It is not clear if the flux goes to zero 
(Ransom, private communication). One possibility is that Pulsar A wind leaks through the interface 
and affects the microphysical process responsible for the generation of radio emission by Pulsar B, 
somewhat similar to the so-called flux transfer events at the day side of Earth magnetosphere (Fahr 
et al. 1986). Such process is prohibited in ideal MHD and should occur through resistive effects 
(e.g. tearing mode). The fact that the magnetized boundary becomes "leaky" and both plasma 
and magnetic field are transported across it, has been amply demonstrated through decades of 
space experiments (e.g. Cowley 1982). The transport occurs either due to microscopic resistive 
instabilities of the surface current (Galeev et al. 1986) or due to dynamic (e.g. RT) instabilities. 

Can Pulsar A affect Pulsar B electrodynamics? We have previously estimated the density of 
Pulsar A wind at Pulsar B to be ua ~ 3 x 10~ 2 A cc. For Pulsar B, the particle flux is N = 
XyJ&ncLB/e 2 = 10 30 As _1 (surface Goldreich-Julian density noj,surf,B = 2.2 x 10 9 cm -3 ), so that 
the density at the interface r m ~ 2 x 10 9 cm is ne ~ 0.27A cc. Thus, if the multiplicity factors are 
equal for both pulsars, the particle density of Pulsar B flow on the open field lines is two orders of 
magnitude larger than the density of Pulsar A wind. Overall, Pulsar A wind can make only a small 
contribution to the particle density in Pulsar B magnetosphere. On the other hand, the assumption 
of equal multiplicity factors may not hold in the gaps of Pulsar B - regions of low density, where 
\b ~ 1. If Pulsar A plasma can get onto field lines passing through Pulsar B gaps, and if A^ > 100, 
this can strongly affect particle acceleration in Pulsar B gaps and radio emission generation. 

Finally, we suggest a possible explanation for sudden increases in brightness of Pulsar B at 
two particular sections in the orbit. The dynamical compression of Pulsar A wind may change 
the condition for radio emission generation if generation takes places at large distances from the 
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pulsar, as has been suggested by Lyutikov et al. (1999a) (see also Kazbegi et al. (1991) For example, 
the growth rate of Cherenkov-drift instability (Lyutikov et al. 1999b) depends sensitively on the 
radius of curvature of magnetic field lines. Due to strong pressure from Pulsar A wind, the open 
field lines of Pulsar B are much stronger curved, than in isolated pulsars, producing larger growth 
rates (growth rates for Cherenkov-drift instability are marginal in isolated pulsars, (Lyutikov et 
al. 1999b)). Cherenkov-drift instability develops at a limited range of radii and produces emission 
which is beamed along the local magnetic field. As Lyutikov et al. (1999a) have argued, emission 
is produced at two locations: in a ring-like region near the magnetic axis and in the region of swept 
back magnetic field lines. In the latter case the radio emission is produced at large angles with 
respect to the magnetic axis (see Fig. 9). We propose that the transient brightening of Pulsar B 
occurs when the line of sight runs parallel to the magnetic field lines in the emission generation 
region located on the bend-back field lines, close to the edge of Pulsar B magnetosphere. 

The immediate implication of the model is that the pulse profile of Pulsar B should change 
with the orbital phase, as the new emission region comes into line of sight at particular parts of the 
orbit. Preliminary data indicate that this is indeed the case (S. Ransom, private communication). 

An important prediction of the Cherenkov-drift instability, which is in stark contrast to the 
the bunching theory of radio emission, is that the emitted waves are polarized perpendicular to the 
plane of the curved magnetic field line. Thus at the swept-back field lines polarization is along the 
axis of rotation. This may be used as a test to distinguish between the two theories. Assuming that 
Pulsar B is an orthogonal rotator with the rotational axis along the normal to the orbital plane, 
the Cherenkov-drift emission will be linearly polarized along the normal to the orbit as well. 

7. Discussion 

In this paper we considered magnetohydro dynamic interaction between relativistic pulsar wind 
and static magnetosphere in double pulsar system PSR J0737-3039. Our main conclusion are: 

• Electromagnetic torque on Pulsar B is increased due to the interaction with Pulsar A wind. 
Depending on the nature of interaction, the magnetic field of Pulsar B can be as low as 
3 x 10 11 G. Still, the geometrical model for the form of the interface predicts the size of the 
interface which is much larger than inferred from the duration of Pulsar A eclipses, unless 
the orbital inclination angle is < 86°. 

• A likely cause of eclipses is synchrotron absorption in the dense shocked Pulsar A wind by low 
energy relativistic electrons. The density of Pulsar A wind should be at the upper allowed 
limit. The presence of such electrons cannot be justified from first principles. 

The major uncertainty of the model is the source of low energy, 7 ~ 30, electrons in the 
sheath. Pulsar A wind cannot be so slow: otherwise induced Compton scattering in the wind will 
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make it unobservable (Sincell & Krolik 1992). One possibility is that such low energy electrons 
are accelerate not at the shock, but at the rotational discontinuity during development of dynamic 
and/or resistive instabilities. Alternatively, mixing of Pulsar A and Pulsar B plasmas, initiated by 
such instabilities, may populate the sheath with weakly relativistic electrons accelerated in Pulsar 
B gaps. If mixing is efficient, particle number density in the sheath may be dominated by Pulsar B 
plasmas. Yet another possibility is that absorption happens on the strongly bend-back open field 
lines of Pulsar B which asymptotically take the form of the sheath (it would be hard to distinguish 
this possibility from absorption in the sheath itself). 

Our calculations are consistent with the possibility that Pulsar A wind experiences a strong 
shock near Pulsar B. This limits the magnetization parameter to a < F 2 (this is a condition that 
the preshock Lorentz factor T is larger than the Alfven speed in the wind). Thus, our results 
do not necessarily imply that Pulsar A wind is kinetically dominated near Pulsar B. (The energy 
required to create a population of low energy electrons is a \ jo fraction of the total wind energy; 
this is also a fraction of energy dissipated in a strongly magnetized shock with u»l, (Kennel & 
Coroniti 1984a).) 

A possible measurement of a may come from high energy observations of the system. Shocks 
in kinetically dominated winds are expected to be strongly dissipative so that a large fraction of the 
incoming energy flux may be radiated in X-rays. The expected luminosity is Lx ~ -L^Afi ~ 3 x 10 30 
erg/sec, where AQ ~ 5 x 10~ 3 is the solid angle of the shocked region seen from Pulsar A. A recent 
detection of a weak X-ray emission (McLaughlin et al. 2004) at a level 2 x 10 30 erg/sec is consistent 
with this scenario. (A confirmation of this result is needed since the count rate was too low and 
the X-ray emission is also consistent with Pulsar A magnetospheric emission.) 

A better understanding of the system should come from the full relativistic MHD modeling of 
the system. The relations for oblique relativistic MHD shocks derived in this paper can be used 
as guidelines and a check to such simulations. Qualitatively, the flow in the sheath is expected 
to become supersonic both due to the changing shock conditions and due to pressure acceleration 
away from the apex point, so that the flow will form a de Laval nozzle. In addition, centrifugal 
forces acting on a flow moving along the curved trajectory (so called Busemann correction) need 
to be taken into account. Another important property of the shocked flow is that even for small 
magnetization of Pulsar A wind a <C 1 magnetic field plays an important role near the contact 
discontinuity and cannot be neglected. Thus, any credible simulation of the interaction must use 
full relativistic MHD formalism. 
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Fig. 1. — Cartoon of the binary pulsar system PSR J0737-3039. The wind of Pulsar A is shocked 
near Pulsar B and forms a sheath. Magnetosphere of Pulsar B is strongly distorted due to inter- 
action with the wind. The system resembles the Earth-Solar wind interaction. On the "day side" 
of Pulsar B magnetosphere (facing Pulsar A) the opening angle of the last closed field line may 
be much larger than in the case of isolated pulsar. In the sheath, the flow near the apex point is 
subsonic, while at larger angles it becomes supersonic both due to changing conditions at oblique 
shocks and due to pressure acceleration along the sheath. The contact surface is also expected to 
be unstable to Kelvin-Helmholtz instabilities. 
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Fig. 2. — Geometry of the interaction. Pulsar A is located at distance D from Pulsar B along the y 
axis. Magnetic moment of Pulsar B may be oriented along any of the axis x, y, z. Pulsar A wind is 
shocked near Pulsar B. At point C on the interface surface, defined by tb{0b), the surface normal 
makes an angle a with re. At the same point, the radial flow from Pulsar A is inclined at angle 4> 
to the surface. The insert shows the eclipsing region (in radians) as viewed from Pulsar A. 
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Fig. 3. — Form of the contact surface in the y — z plane for different orientations of Pulsar B 
magnetic dipole. Pulsar B is located at the origin, Pulsar A is located at y = 1. Upper curves for 
magnetic field of Pulsar B assuming standard spin down luminosity (Bb = 1-5 x 10 12 G). Solid 
lines are for \i = fi x (in the plane of orbit along the Pulsar B orbital velocity), dashed lines are 
for fjb = n y (along the line connecting the two pulsars) and dotted lines for n = ji z ( perpendicular 
to the orbital plane; observationally, large fj, z is excluded since we do see the pulses from Pulsar 
B). Asymptotic opening angle is ~ 0.11 — 0.13 rad. Lower solid curves for magnetic field of Pulsar 
B with modified torques Bb = 4.7 x 10 11 (Eq. 10) (asymptotic opening angle is ~ 0.09 rad) and 
Bb = 3 x 10 11 G (Eq. 12) (asymptotic opening angle is ~ 0.07 rad). Long dashed line is the line 
of sight (inclination 87.7°). Dash-dot line is the light cylinder radius (it looks non spherical due to 
different axial scales). 
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Fig. 4. — Angle of attack <p as a function of polar angle. Negative 9b refer to parts of the interface 
facing Pulsar A. Labeling of the curves is the same as in Fig. (3). 
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Fig. 5.- 



Geometry of oblique shock flow. 




Fig. 6. — Dependence of the deflection angle % on the angle of attack (f> for initially cold (pi = 0) 
fluid for different u\ and T = 4/3. u\ — > corresponds to non-relativistic case. 



-23- 




Fig. 7. — Dependence of the deflection angle % on the angle of attack 4> f° r ultra-relativistic 
MHD shocks ui > 1 and different magnetization parameters a. Dotted lines gives the maximum 
deflection angle (Eq. (33)), dashed line demarcates the subsonic shocked flow (to the right) and 
supersonic flow (Eq. (34)). For fluid shocks the maximum deflection angle Xmax = 30° is reached 
at 4>max = 60°. This is also the point when the shocked flow becomes subsonic. 
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Fig. 8. — Post-shock Mach number as a function of polar angle. (a) hydrodynamic flow 
(7 = (post-shock four- velocity is M2/V2), (b) magnetized wind a = 1 (post-shock four- velocity 
■v/(l + 3<r)/2M 2 ). Labeling of the curves is the same as in Fig. (3). 
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Fig. 9. — (a) Fig. 5 from Lyutikov et al. (1999a) showing the location of radio emission generation 
regions in isolated pulsars (field structure is that of a vacuum rotating dipole). "The location 
of the Cherenkov-drift emission depends sensitively on the curvature of magnetic field line. Two 
possible locations of the Cherenkov-drift emission are shown: ringlike near the magnetic axis and 
in the region of swept field lines." (b) Possible modification due to interaction with Pulsar A 
wind. The dim component of emission is not affected by the Pulsar A wind and is produced at all 
phases of Pulsar B rotation, while the bright component is produced only at the side of Pulsar B 
magnetosphere facing Pulsar A. 



